import sympy
# 直接求
x,y,z=sympy.symbols('x y z')
fx=1
fy=1
fz=x
res=sympy.integrate(sympy.integrate(sympy.integrate(fz,(z,0,1-x-2*y))*fy,(y,0,(1-x)/2))*fx,(x,0,1))
sympy.pprint(res)